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摘要 利用 廊坊 站 (40.0" N, 116.3° E) 钠 层 测 温 测 风 激 光 雷 达 2011 年 至 2013 年 共 约 82h 的 钠 
原子 数 密度 和 垂直 风 观 测 数 据 , 首次 分 析 了 廊坊 地 区 中 间 层 顶 区 域 大 气 重力 波 耗 散 引 起 的 钠 
原子 输送 作用 . 分 析 得 到 ，90~100 km 重力 波 耗 散 引 起 的 平均 钠 原 子 垂直 通 量 整体 为 负 ， 钠 
原子 向 下 输送 ， 在 93 km 处 出 现 最 大 负 值 为 -1.47x10* m”m/s，85~90 km 平均 钠 原 子 垂直 通 
> 量 为 正 ， 钠 原子 向 上 输送 , 但 通 量 值 随 高 度 递减 . 钠 原 子 垂 直通 量 方向 在 90 km 处 发 生 转变 . 
垂直 通 量 随 高 度 的 变化 造成 钠 原 子 汇聚 ， 汇 聚 效应 引起 的 平均 钠 原 子 产生 率 最 大 值 在 91km 
处 达到 了 1.4x108m-3/h， 该 值 超过 了 相同 高 度 上 理论 模拟 的 流星 烧 蚀 注入 引起 的 钠 原 子 产 生 
率 峰 值 ， 说 明 它 对 钠 层 结构 的 形成 具有 重要 的 贡献 . 分 析 结 果 与 美国 的 观测 结果 相 比 ， 平 均 
钠 原子 产生 率 峰 值 大 小 相近 , 但 出 现 高 度 不 同 , 说 明 大 气 重力 波 耗 散 引起 的 物质 输送 具有 较 
大 的 地 域 变化 特征 . 研究 结果 可 为 大 气 物质 输送 理论 的 完善 以 及 大 气 金属 层 物理 模式 的 改 
进 提供 观测 事实 参考 
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中 间 层 项 区 域 是 大 气 重 力 波 活 动 和 耗 散 非常 显著 的 区 域 (Gardner et al., 2002; Zhao Y et 
al., 2003). 大 气 重 力 波 引 起 大 气 粒 子 在 局 部 平衡 位 置 附 近 发 生 振荡 ， 如 果 大 气 重力 波 发 生 
耗 散 ， 振 荡 的 大 气 粒 子 不 能 回 到 其 平衡 位 置 ， 从 而 产生 对 粒子 的 输送 作用 (Walterscheid and 
Hocking, 1991; Walterscheid, 2001); 同时 ， 如 果 粒 子 离开 平衡 位 置 后 与 其 他 物质 发 生化 学 
反应 ， 也 会 产生 粒子 输送 作用 ( 徐 寄 遥 等 ，2001)， 这 种 输送 作用 是 大 气 波动 与 化 学 过 程 耦合 
引起 的 (Walterscheid and Schubert, 1989). 大 气 物 质 的 垂直 输送 是 中 间 层 顶 大 气 结构 和 动力 
学 研究 的 前 沿 课题 . 目前 ， 在 全 球 大 气 层 模式 (如 WACCM) 研 究 中 ， 仅 考虑 滑 流 混合 机 制 的 
物质 垂直 输送 ， 没 有 考虑 大 气 重 力 波 耗 散 引起 的 垂直 输送 机 制 (Liu and Gardner，2004)， 模 
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式 对 中 间 层 顶 钠 层 的 模拟 结果 与 观测 结果 还 有 较 大 的 差别 (Plane，2003). 大 气 重 力 波 耗 散 现 
象 对 大 气 各 种 物质 成 分 垂直 输送 过 程 中 的 贡献 , 是 一 个 知之 不 多 的 、 竺 解决 和 参数 化 的 科学 


问题 . 

钠 层 测 风 测 温 激光 雷达 提供 了 从 观测 手段 上 揭示 大 气 重 力 波 耗 散 引起 物质 垂直 输送 作 
及 其 影响 的 方式 . 该 激光 雷达 可 同时 探测 中 间 层 顶 区 域 (80~110 km) 的 大 气 三 维 风 场 、 温 度 
及 钠 原 子 数 密 度 . 中 间 层 顶 钠 原 子 因 具 有 较 长 的 寿命 和 相对 稳定 的 状态 可 作为 分 析 大 和 气 重 
力 波 耗 散 输送 作用 的 示 踪 物 (Bowman et al., 1969). Liu 和 Gardner(2004) 首 次 利用 钠 层 测 风 测 
温 激光 雷达 在 SOR 站 (35.0"N，106.5°W) 的 观测 数据 估算 了 大 气 重力 波 耗 散 对 钠 原 子 输送 的 
显著 效应 ， 并 与 涡 旋 输送 结果 进行 对 比 ， 结 果 显示 两 者 对 钠 原 子 输送 贡献 相当 . Liu 和 
Gardner(2005) 对 SOR 和 Maui 站 (20.7*N，156.3°W) 钠 激光 雷达 观测 数据 进行 对 比 研究 ， 发 
现 重 力 波 耗 散 引 起 的 输送 作用 随地 理 位 置 而 变化 . Gardner 和 Liu(2010) 利 用 SOR 的 长 期 观测 
数据 ， 得 到 了 重力 波 引 起 的 钠 原 子 输送 作用 及 多 种 输送 速度 的 季节 变化 规律 
目前 ， 虽 然 上 述 研究 结果 对 中 间 层 顶 重力 波 -大 气 成 分 耦合 过 程 的 理解 帮助 很 大 ， 但 离 
我 们 全 面 掌握 大 气 重 力 波 耗 散 输 送 作用 、 全 球 变化 特性 、 参 数 化 方案 等 还 有 很 远 的 距离 ， 需 
要 更 多 台 站 的 观测 分 析 结 果 来 支撑 . 本 文 首次 给 出 利用 位 于 东 半 球 的 中 科 院 廊坊 临近 空间 
环境 野外 科学 观测 研究 站 (40.0°? N, 116.3°E) 钠 层 测 风 测 温 激 光 雷 达观 测 数据 分 析 大 气 重力 波 
耗 散 引起 的 钠 原 子 输送 作用 的 结果 ， 并 与 美国 SOR 与 Maui 的 观测 结果 进行 比较 
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2 观测 设备 和 数据 来 源 


廊坊 站 钠 层 测 风 测 温 激光 雷达 是 中 国 科学 院 国 家 空间 科学 中 心 自主 研发 的 先进 探测 设 


备 (Yan et al., 2009, Hu etal., 2011; 徐 丽 等 , 2010), 其 探测 原理 与 美国 的 同类 激光 雷达 (Bills 


etal., 1991; She and Yu, 1994; Gardner and Yang，1998) 相 同 . 该 激光 雷达 发 射 三 束 激光 ， 
分 别 指向 天 顶 方向 、 天 项 偏 东 20*、 天 顶 偏 北 20"， 采 用 3 个 望远镜 分 别 接收 以 上 三 个 方向 
的 大 气 回 波 信号 ， 以 同时 探测 80~105 km 大 气 三 维 风 场 、 温 度 和 钠 原 子 数 密度 . 系统 参数 见 
#1. 

该 激光 雷达 自 2010 年 底 以 来 开始 常规 观测 ， 其 温度 、 风 场 和 钠 原 子 密度 等 观测 数据 已 
于 中 间 层 顶 温度 结构 (Hu et al.，2011)、 准 单 色 大 气 重力 波 活 动 特性 ( 王 博 等 ，2016)、 突 发 
内 层 观测 分 析 ( 程 永 强 等 ,2016) 以 及 大 气 模型 验证 等 研究 工作 . 本 文选 取 2010~2013 年 共 82 
h 的 夜间 垂直 风 场 和 钠 原 子 数 密度 观测 数据 进行 大 气 重 力 波 耗 散 引起 的 钠 原 子 垂直 输送 特 
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征 统计 分 析 ， 观 测 数 据 的 时 间 分 辨 率 15 min, FEATHER 1 km， 最 长 一 天 的 观测 时 长 h, 
平均 观测 时 长 为 4.5 h. 


表 1 钠 菊 光 多 普 惑 激光 雷达 系统 参数 


Table 1 System parameters of sodium fluorescence Doppler lidar 


名 称 数值 
激光 波长 589.158 nm 
调制 频率 +630 MHz 

激光 平均 功率 ~1W 

激光 重复 频率 30 Hz 

激光 脉冲 长 度 ~7 ns 
激光 发 散 角 ~1 mrad 
望远镜 口径 Im 


3 分 析 方 法 


为 了 从 观测 数据 中 统计 得 到 大 气 重 力 波 耗 散 引起 的 钠 原 子 垂直 输送 , 首先 需要 提取 大 气 
重力 波 引起 的 垂直 风 扰动 、 钠 原子 数 密度 扰动 , 然后 用 大 气 重力 波 扰动 量 统计 各 高 度 的 上 平 
均 钠 原子 垂直 通 量 、 平 均 钠 原子 产生 率 及 其 不 确定 度 . 下 面 介绍 相关 分 析 方 法 . 


3.1 大 气 重力 波 扰 动量 提取 


利用 每 晚 观 每 个 高 度 上 的 垂直 风 和 钠 原 子 数 密度 观测 数据 , 组 成 时 间 序 列 , 原始 数据 减 
去 该 时 间 序列 的 线性 趋势 ,得 到 的 残 差 值 作为 垂直 风 扰动 值 和 钠 原 子 数 密度 的 扰动 值 . 由 于 
所 采用 的 观测 数据 垂直 分 状 率 为 1 km、 时 间 分 状 率 为 15 min， 平 均 观测 时 间 4.5 h， 所 以 垂 
直 风 和 钠 原 子 数 密度 扰动 值 包含 了 平均 周期 为 30 minh 以 内 、 垂 直 波 长 为 2km 以 上 的 大 
气 重 力 波 信息 . 所 有 的 数据 首先 经 过 三 倍 标准 偏差 法 检验 以 剔除 毛刺 点 . 图 1 为 2010 年 12 
月 28 日 大 气 重力 波 重 直 风 扰 动 提取 的 示例 . 
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图 1 大 气 重 力 波 扰动 提取 


Fig. 1 Extraction of atmospheric gravity wave perturbations 


3.2 大 和 气 重力 波 扰动 数据 验证 


为 了 验证 提取 的 大 气 重 力 波 扰动 值 的 可 靠 性 ， 需 要 检验 扰动 值 的 分 布 情况 . KAEN 
活动 可 以 认为 是 大 气 中 的 一 个 随机 现象 , 其 引起 的 地 球 物理 量 的 变化 具有 随机 分 布 特性 . 
中 心 极限 定理 (Central Limit Theorem) 可 知 ， 随 机 变量 分 布 渐 近 于 高 斯 分 布 ， 大量 的 随机 变 
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积累 分 布 会 逐 点 收敛 到 高 斯 分 布 . 根据 观测 数据 提取 的 垂直 风 和 钠 原子 数 密度 扰动 值 是 


实 值 和 测量 误差 的 和 , 测量 误差 的 来 源 是 光子 噪声 ， 虽 然 光 子 噪声 属于 泊 松 分 布 , 但 泊 松 
布 在 大 量 粒子 数 时 趋向 于 正 态 分 布 ， 因 此 ， 激 光 雷 达 的 大 气 重力 波 扰动 值 ( 含 测量 误差 ) 符 
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a 
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高 斯 随机 分 布 过 程 . 在 实际 大 气 中 ， 当 相互 独立 的 波动 的 数量 足够 大 的 时 候 ， 垂 直 风 扰动 、 


钠 原 子 数 密度 扰动 及 测量 误差 就 应 该 符合 均值 为 0 的 高 斯 随机 分 布 模型 . 
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图 2 (a) 钠 原子 数 密度 扰动 分 布 及 高 斯 拟 合 
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，(b) 垂 直 风 扰动 分 布 及 高 斯 拟 合 


Fig. 2 (a) Distribution of Na density perturbations and Gaussian fitting; (b) Distribution of 


vertical wind perturbations and Gaussian fitting. 


对 提取 到 的 垂直 风 扰 动 和 钠 ) 


原子 数 密度 扰动 按照 分 布 范围 


钠 原 子 数 密度 扰动 和 垂直 风 扰 动 均 属 于 正 态 
垂直 风 扰 动 主要 分 布 在 土 10m/s 的 范 
结果 较 好 , 原因 是 大 气 波动 是 垂直 风 
响 钠 原子 数 密度 变化 的 因素 要 比 垂直 风 复 杂 的 多 ， 除 受 波动 影响 外 ， 还 会 
JEH. 总 体 上 说 ， 二 者 的 扰动 值 分 布 都 较 好 的 符合 高 其 
表明 钠 原 子 数 密度 和 垂直 风 扰 动 数据 是 可 靠 的 . 


合 结果 比 钠 原 子 扰动 的 拟 合 


化 学 反应 、 中 性 


分 布 模型 ， 


化 学 反应 等 因素 共同 


3.3 大 气 重力 波 耗 散 引 起 的 钠 原 子 垂直 通 量 


论 上 ， 钠 原子 的 垂直 通 


WPw 表示 平均 背景 垂直 风 引 起 的 钠 ) 


均 垂 直通 量 ， 


垂直 通 量 w'y 是 本 文 研究 的 对 象 . 值得 注意 的 是 非 耗 散 重 力 波 产生 的 w 和 AD 


关系 为 正 交 ，w' py, 为 零 ， 所 以 非 耗 散 重 力 波 不 会 引起 垂直 通 量 ; 


散 时 ，W' 和 p',, 原本 正 交 的 相位 关系 被 改变 ，w' Pw DEF, 


引起 垂直 通 量 . 


w' 和 p' 分 别 表示 垂直 风 扰动 和 钠 原 


原子 垂直 通 


量 是 垂直 风 和 钠 原 子 数 密度 的 乘积 ， 如 公式 (1)， 


进行 高 斯 拟 合 ， 如 图 


2 所 示 . 


分布, 钠 原子 数 密度 扰动 主要 集中 在 士 3x10"m ， 
BAY. 从 图 2 中 的 分 布 结构 可 知 ， 垂 直 风 扰动 的 高 斯 拟 


产生 扰动 的 动力 来 源 , 而 影 
受到 流星 注入 、 光 


右边 第 一 项 


量 ， 第 二 项 w'pw 表 示 大 气 波动 引起 的 平 


原子 数 密度 扰动 . 大 气 重力 波 引 起 的 平均 


其 相位 


而 当 大 气 重力 波 发 生 耗 


所 以 大 气 重力 波 耗 散 能 够 


上 节 中 从 观测 数据 中 提取 得 到 的 大 气 重 力 波 扰动 包含 测量 误差 ,所 以 利用 观测 数据 统计 


大 气 重 力 波 引起 的 钠 原 子 垂直 通 量 如 公式 (2) 所 示 ，CA 表示 重力 波 引起 的 钠 原 子 垂直 通 量 ， 


Aw FI Apya 表示 垂直 风 和 钠 原 子 数 密度 的 测量 误差 . 公式 (3) 是 大 气 重 力 波 引 起 的 平均 钠 原 


子 垂直 通 量 . 每 晚 观 测 的 大 气 重力 波 扰动 和 测量 误差 的 平均 值 为 0， 即 ww'、PWw 、Aw 和 


Apy 均 为 零 ， 因 此 由 公式 G3) 可 知 钠 原 子 垂 直通 量 统计 结果 不 受 激光 雷达 测量 误差 的 影响 . 
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对 于 激光 雷达 每 一 个 晚上 的 连续 观测 数据 , 平均 钠 原 子 垂直 通 量 的 方差 估计 式 如 公式 (4) 


所 示 ， 其 中 7 是 观测 时 长 ，L=1 km 是 垂直 高 度 分 辨 率 ，Azor =lkm 为 钠 原 子 垂直 通 量 相 关 


AZ, Ator =15 min 为 钠 原子 垂直 通 量 相关 时 间 ( 参 量 的 取 值 及 方差 推导 过 程 请 详 见 Gardner 


and Yang，1998)， 观 测 到 的 垂直 风 扰 动 平方 的 平均 值 就 表示 为 (ww ) =(w') +(Amw) . 
根据 每 晚 提取 的 垂直 风 及 钠 密度 扰动 数据 ， 利 用 公式 (3) 可 统计 得 到 每 晚 各 个 高 度 上 大 

© 气 重力 波 引起 的 钠 原 子 垂直 通 量 . 根据 每 晚上 的 观测 时 间 长 度 z ,利用 公式 (4) 可 统计 得 到 每 

晚 各 个 高 度 上 平均 钠 原 子 垂直 通 量 的 方差 . 用 公式 (5) 计 算得 到 所 有 晚上 各 个 高 度 上 的 平均 


钠 原 子 垂直 通 量 方差， 对 公式 (53) 的 统计 结果 进行 开 方 计算 ， 得 到 平均 钠 原 子 垂直 通 量 的 标 


准 差 作为 各 个 高 度 上 平均 钠 原子 垂直 通 量 的 不 确定 度 参数 . 


W: Pha =(w+w')(Pya + Pya ) 
=WPya + W' Pya ‘+ Ppa W E W' Pya (1) 
=WPna +W' Pya "+ PygW'+W' Pya 


= WPywa +W' Pya ' 
GF =(w'+Aw):(p'y,+Apy,) (2) 


GF =(w'+ Aw) (p'u + APy,) 
=W'P' yat P'ya AWt+w'Ap,, + AWAPya 


{°° (3) 
=W'P'va + P ng AWtw'Ap,, + AwAPy, 
= W'P'ya 
r/2 LI2 
Var(GF)=— | f [oo+rAwom+Aow-wPwPadidc 
T —t/2-L/2 (4) 


heer or (oF (oy | stp] 


L 


_ N = 
Var(GF) ,ww = > Var(GF), (5) 
i=l 


3.4 垂直 通 量 引起 的 钠 原 子 产 生 率 


中 间 层 顶 区 域 大 气 重 力 波 耗 散会 引起 钠 原 子 垂直 通 量 随 高 度 变 化 , 钠 原 子 垂 直通 


量 的 梯 


度 会 导致 不 同 高 度 的 钠 原 子 汇聚 或 者 流失 .由 重力 波 耗 散 引起 的 钠 原子 数 密度 增加 


或 减 小 


可 以 用 本 地 钠 原 子 产 生 率 (Na Production Rate) 来 表示 ， 如 公式 (6)， 其 中 已 ,表示 钠 原 


R, Pu 为 正 表示 本 地 钠 原 子 受 重力 波 耗 散 作用 发 生 堆积 ， 数 密度 增加 ， Py, 为 负 即 


率 ， 表 示 本 地 钠 原 子 流失 ， 数 密度 减 小 (Gardner and Liu, 2010). 公式 (7 中 Var(P,) 


晚 各 个 高 度 上 的 平均 原子 产生 率 的 方差 , T 是 观测 时 长 ,L 是 垂直 高 度 范围 , 厂 为 标高 


为 钠 原 子 垂直 通 量 相关 高 度 ，Atos 为 钠 原 子 垂直 通 量 相关 时 间 . 


_ OWw'P'wa 


子 产生 


pe — 和 一 


ARR EF 


» Azor 


P = (0) 


NM 
a Oz 


(P'wa) +(APna) 


mel 


Var(P,, 


a 


2 (wY +(Aw)? | 
=a) | al dL 
L L T H 


(7) 


— 1 Ww — 
Var(Py, ieee = N? > Var(P。 ); (8) 
i=l 


根据 公式 (3) 得 到 每 晚 平 均 钠 原子 垂直 通 量 的 统计 结果 ， 利 用 公式 (6) 统 计 得 到 每 晚 各 个 
高 度 上 的 平均 钠 原 子 产生 率 . 利用 公式 (7) 可 统计 得 到 每 晚 各 个 高 度 上 的 平均 钠 原 子 产 生 率 


的 方差 . 利用 公式 (8) 统 计 得 到 所 有 晚上 各 个 高 度 上 的 平均 钠 原 子 产生 率 的 方差 ， 对 


公式 (8) 


的 统计 结果 进行 开 方 计算 , 得 到 各 个 高 度 上 平均 钠 原 子 产 生 率 的 标准 差 作为 各 个 高 度 上 平均 


钠 原 子 产生 率 的 不 确定 度 参数 . 


3.5 测量 误差 的 影响 


钠 激 光 雷 达 测 量 的 垂直 风速 度 和 钠 原 子 数 密度 数据 都 含有 测量 误差 . 测量 误差 对 估计 
值 的 影响 大 小 与 平均 钠 原 子 垂 直通 量 和 产生 率 的 统计 结果 是 否 可 靠 十 分 相关 . 从 平均 钠 原 
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子 垂直 通 量 估计 公式 (3) 和 平均 钠 原 子 产 生 率 估计 公式 (6) 来 看 ， 测 量 误差 不 影响 这 些 估 计 值 . 
所 以 利用 上 述 方法 估计 的 钠 原 子平 均 垂直 通 量 和 平均 产生 率 是 可 靠 的 . 
从 公式 (0) 和 (7) 来 看 ， 每 晚 统计 的 平均 钠 原 子 垂直 通 量 和 钠 原 子 产生 率 的 方差 不 仅 受 测 


由 


量 误差 Aw Fl Apy, 的 影响 ， 而 且 受 大 气 重 力 波 扰动 WwW 和 ww 发 生地 球 物理 学 变化 的 影响 


(Fritts, 2000). 为 减 小 二 者 的 统计 方差 , 使 得 平均 垂直 通 量 和 平均 产生 率 有 确定 的 物理 意义 ， 
我 们 采取 增加 观测 数据 时 间 尺 度 和 空间 尺度 的 办 法 . 本 文 我 们 采用 了 82h 的 时 间 尺 度数 据 和 
垂直 Ikm 的 空间 尺度 数据 ， 给 出 有 统计 意义 的 结果 . 


4 结果 和 讨论 


4.1 钠 层 结构 


105 T T T T T T 


85+ 4 


80 | | | | 1 
0 0.5 1 1.5 2 2.5 3 3.5 


Na Density m°’) x10 
图 3 RAF AAE T A RE 
Fig. 3 Mean Na atoms density profile at LangFang 


KAET AAE PAR 8 PS ai a PEP TD, WE 3. 廊坊 地 
平均 钠 原 子 数 密度 垂直 剖面 呈 准 高 斯 分 布 结 构 ,， 钠 层 峰 值 高 度 为 90km, 质心 高 度 为 91km， 
平均 数 密 度 最 大 值 为 3.17x10”m”. 中 间 层 顶 金 属 钠 的 来 源 是 流星 在 大 气 中 烧 蚀 ， 此 外 钠 原 
子 数 密度 变化 还 受 电离 层 、 光 化 学 反应 及 大 气 波动 因素 的 共同 作用 (Hedin and Gumbel, 2011). 
太阳 辐射 对 中 间 层 项 上 部 钠 离 子 参与 的 光化学 反应 产生 直接 影响 (Mcneil et al.，2002); FA] 


xl 
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层 顶 活性 物质 如 O3. OH 的 浓度 变化 影响 钠 原 子 和 及 其 化 合 物 所 参与 的 中 性 化 学 反应 (Plane 
and Helmer，1994)， 对 钠 层 中 下 部 的 钠 原 子 和 钠 化 合 物 影响 明显 ; 重力 波 、 潮 汐 等 大 气 波动 
和 小 尺度 的 涡 旋 对 钠 原 子 运动 和 输送 提供 了 能 量 ， 造 成 钠 原 子 垂直 输送 和 混合 (Fritts and 


Alexander, 2003). 


这 些 因素 相互 联系 ， 相 互 作用 共同 作用 形成 准 高 斯 分 布 的 钠 层 结构 形态 


4.2 大 和 气 重力 波 耗 散 引起 的 钠 原 子 垂直 输送 
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Altitude (km) 


8 


—— Vertical Na Flux 


i 
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 
Vertical Na Flux (m? m/s) x10 


图 4 廊坊 站 平均 钠 原 子 垂直 通 量 廓 线 及 其 标准 差 


Fig.4 Mean profile of vertical Na flux and its standard deviation at LangFang 


图 4 Ah SPA eae EAR Re, RAR ZAI 5km WA AMET PIA, E 


线 正 值 表示 垂直 
通 量 的 标准 差 表示 


向 上 运动 ， 负 值 表示 垂直 向 下 运动 ， 误 差 棒 为 不 确定 度 ， 由 平均 钠 原 子 垂直 
示 ， 误 差 棒 的 大 小 不 影响 钠 通 量 最 大 值 和 最 小 值 正 负 的 判断 . 图 中 明显 看 到 


钠 通 量 在 在 90 km 钠 层 峰值 高 度 ( 见 图 3) 附 近 发 生 方向 转变 . 90 km 以 下 ， 钠 通 量 为 正 ( 向 上 
输送 ); 90 km UE, 钠 通 量 为 负 ( 向 下 输送 ), 在 93 km 处 达到 最 小 值 -1.47x108m3 m/s; 97 km 


附近 ， 钠 通 量 接 


近 于 0. 
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图 5 廊坊 站 钠 原 子 产 生 率 廓 线 及 其 误差 


Fig. 6 Mean profile of Na production rate and its standard deviation at LangFang 

图 5 给 出 了 钠 原 子 产 生 率 剖 面 及 其 误差 棒 . 产生 率 廓 线 经 过 5 km 汉 明 窗 进行 平滑 . 误差 
棒 为 不 确定 度 ,， 由 平均 钠 原子 产生 率 标准 差 表示 ,误差 棒 的 大 小 不 影响 产生 率 最 大 值 和 最 小 
值 正 负 的 判断 . 在 93.5 km 以 下 ， 钠 原子 产生 率 为 正 值 ， 说 明 此 高 度 范围 内 重力 波 耗 散 作用 
使 得 钠 原 子 数 目 在 单位 时 间 内 产生 净 增 长 ， 平 均 钠 原子 产生 率 最 大 值 在 91 km 处 达到 了 
1.47x10* mh. 钠 层 的 质心 高 度 与 钠 原 子 产生 率 最 大 值 高 度 一 致 ， 这 可 以 说 明 大 气 重 力 波 耗 
散 输 送 作 用 使 得 钠 原子 在 91 km 质心 高 度 附 近 产 生 堆 积 , 大 气 重力 波 耗 散 引起 的 钠 原 子 输送 
作用 可 能 对 钠 层 质心 高 度 附 近 的 峰值 结构 特征 的 形成 起 着 重要 的 作用 . 


K 3 给 出 了 廊坊 站 (40.0"N，116.3°E) 与 美国 SOR(35.0?N，106.5°*W) 和 Maui(20.7°N， 


156.3"W) 观 测 站 (Liu and Gardner，2005) 的 观测 统计 结果 的 比较 . 表 3 显示 三 个 地 区 钠 原 子 
垂直 通 量 最 小 值 和 钠 原 子 产生 率 最 大 值 基本 相当 ， 但 是 其 平均 扩 线 趋势 存在 差异 ，SOR 地 
区 和 Maui 地 区 在 85~100 km 的 垂直 钠 通 量 基 本 为 负 ， 最 大 负 值 分 别 出 现 在 87 km 和 88 km 
附近 ， 钠 原子 产生 率 最 大 值 在 86 km 附近 ; 廊坊 地 区 85~90 km 钠 通 量 为 正 向 ， 最 大 负 值 出 
现在 93 km 附近 , 钠 原 子 产 生 率 最 大 值 在 91 km 附近 . 三 个 观测 点 统计 结果 存在 不 同 的 原因 
可 能 是 由 于 地 理 位 置 等 条 件 差 异 造成 ， 目 前 全 球 对 于 钠 原子 垂直 输送 的 观测 还 是 较 少 ,其 中 
的 差异 问题 还 需 更 多 数据 分 析 来 了 解 . 
表 3 廊坊 、SOR 和 Maui 三 站 钠 原 子 垂直 通 量 和 钠 原 子 产 生 率 比 较 


Table 3 Comparison with vertical Na flux and Na production rate between Langfang, 


SOR and Maui stations 
观测 点 垂直 通 量 最 大 负 值 及 高 度 产生 率 最 大 负 值 及 高 度 
廊坊 
-1.5x10°m* m/s ~93 km 1.4x10s m3/h ~91km 
(40.0°N, 116.3°E) 
SOR 
-2.6x10°m* m/s ~87 km 2.8x10°m°/h ~86 km 
(35.0°N, 106.5°W) 
Maui 
-0.8x10°m* m/s ~88 km 1x108m3h ~86km 


(20.7°N, 156.3°W) 


4.3 讨论 


一 般 认 为 钠 层 准 高 斯 分 布 的 结构 主要 受 五 部 分 因素 影响 , 第 一 部 分 是 流星 烧 蚀 注入 效应 ， 

由 于 钠 层 的 高 度 与 流星 烧 蚀 高 度 相 吻合 ,所 以 流星 伐 蚀 产生 的 钠 原 子 可 直接 注入 90 km 附近 
的 钠 层 峰值 (Plane，2004); 第 二 部 分 是 离子 化 学 过 程 ，90 km 以 上 钠 离子 化 学 过 程 占 主导 地 
位 ， 部 分 Na 与 NO 和 0; 发 生 一 系列 离子 反应 从 而 变 为 Na ， 部 分 Na 直接 受 太 阳 辐 射电 离 
成 为 Na*+，Na! 数 密度 峰值 在 100 km 附近 (Mcneil et al., 2002); 第 三 部 分 是 中 性 化 学 过 程 ， 
85~90 km 以 下 中 性 化 学 反应 占 主导 地 位 ，Na 原子 与 03、H2，CO; 发 生 系列 反应 变 为 状态 稳 
定 的 NaHCO3， 并 向 低层 沉积 (Plane，2003)，NaHCO; 数 密度 峰值 在 85km 附近 . 根据 模式 
计算 得 到 Na* 和 NaHCO3 的 数 密度 是 Na 原子 数 密度 峰值 0.2~0.3 倍 , 离子 化 学 和 中 性 化 学 过 
程 对 钠 层 峰值 的 页 献 有 限 (Plane，2004); 第 四 部 分 是 涡 旋 输送 过 程 ，85~95 km 涡 旋 起 均匀 
混合 作用 , 混合 作用 造成 钠 原子 向 低层 大 气 输送 ,对 钠 层 准 高 斯 及 峰值 结构 的 形成 起 负 作用 
(Liu and Gardner, 2004). 第 五 部 分 是 大 气 重力 波 输送 过 程 , 大 气 重 力 波 耗 散 引起 的 钠 原 子 垂 
通 量 会 造成 钠 原 子 在 垂直 方向 上 产生 堆积 或 流失 , HE 92km 处 重力 波 耗 散 引起 的 钠 原 子 
产生 率 与 流星 烧 蚀 注入 引起 的 钠 原 子 产生 率 二 者 相当 (Gardner and Liu，2010)， 由 此 可 见 大 
气 重 力 波 耗 散 引起 的 钠 原 子 扒 积 效应 与 流星 烧 蚀 产生 的 钠 原 子 注入 效应 对 钠 层 峰 值 结 构 形 
成 具有 重要 的 作用 . 
本 文 研究 结果 显示 廊坊 地 区 大 气 重力 波 耗 散 引 起 的 输送 作用 造成 钠 原 子 在 质心 高 度 91 


km 附近 汇聚 . Plane(2004) 利 用 NAMOD 模式 和 LDEF(Long Duration Exposure Facility) 观 测 数 


上 


据 ， 分 析 了 流星 在 MLT 区 域 的 烧 蚀 高 度 、 
平均 速度 为 18 km/s 的 流星 在 92 km 烧 蚀 注入 引起 的 钠 原子 产生 率 峰 值 为 2.5x10 mh. 廊坊 
91km 的 峰值 产生 率 达 到 了 1.4x10sm>h， 超 过 流星 烧 蚀 注入 引起 的 产生 率 
峰值 . 上 述 观 测 结果 进一步 说 明 大 气 重 力 波 耗 散 对 钠 层 的 分 布 及 峰值 形成 的 贡献 作用 不 可 


站 统计 结果 显示 


忽视 ， 当 前 没有 考虑 这 种 输送 机 


本 文 首 次 采 


中 科 院 国家 空间 科学 中 ， 


和 垂直 风 观 测 数据 , 分 析 了 中 间 层 项 重力 波 


烧 蚀 效率 及 烧 蚀 产生 的 钠 原 子 注 入 量 ， 结 果 显 示 


判 的 大 气 模型 需要 完善 . 


心 序 坊 站 钠 层 测 温 测 风 激光 雷达 的 钠 原 子 数 密度 
耗 散 引起 的 钠 原 子 垂直 输送 作用 及 其 对 钠 层 结构 


的 影响 , 廊坊 地 区 大 气 重力 波 耗 散 引 起 的 钠 原 子 垂直 通 量 在 90 km 以 下 为 正 ， 钠 原子 向 上 和 输 
送 ，90 km 以 上 为 负 ， 钠 原子 向 下 输送 ，93 km 处 达到 最 大 负 值 -1.47x10; m? m/s; 91 km 钠 


层 质 心 处 的 钠 原 


子 产 生 率 达 到 最 大 值 1.47x108 m3h， 大 于 理论 计算 的 流星 烧 蚀 注入 引起 的 
钠 原子 产生 率 峰 值 ， 说 明 大 气 重 力 波 耗 散 引起 的 输送 作用 对 整个 钠 层 结构 形成 有 重要 影响 . 
廊坊 的 观测 结果 与 美国 SOR 和 Maui 两 地 观测 结果 相 比 ， 基 本 一 致 ， 但 是 钠 原 子 垂直 通 量 


方向 转换 高 度 和 产生 率 峰 值 高 度 不 同 ， 说 明了 大 气 重 力 波 耗 散 输送 作用 具有 地 域 性 特征 . 
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Abstract A total of about 82 h observation data derived from a sodium fluorescence Doppler lidar 
of National Space Science Center (NSSC) at LangFang station, Hebei province, were used to 
study the transport of Na atoms induced by gravity wave dissipation in mesopause region. The 
mean vertical Na flux was mostly negative between 90 km and 100 km and the maximum was 
-1.47x10° m? m/s at 93 km, Na atoms were transported downward, the mean vertical Na flux was 
positive from 85 km to 90 km and Na atoms were transported upward. The mean Na production 
rate was positive from 85 km to 93.5 km and the maximum was 1.4x10° m°/h at 91 km, from 93.5 
km to 100 km, the mean Na production rate was mostly negative. Gravity wave dissipation make 
Na atoms gathered at near 91 km and make the mean Na production rate significantly increase 
between 88 km and 92 km. The mean Na production rate induced by gravity wave dissipation 
significantly exceeds the Na production rate induced by meteoric input at near 91 km. It indicates 
gravity wave dissipations have an important contribution to the formation of the Na layer structure. 
Compared with the SOR and Maui observations, magnitude of Na flux and Na production are 
comparable in three stations, but the mean profiles of Na flux and Na production rate in three 
stations are not quite the same, it indicates geographical factors make differences to atmospheric 
material transport in different places. These results can provide guidance to parameterization and 


evolution of atmospheric metal model in mesopause region. 


Keywords: Sodium fluorescence Doppler lidar; Gravity wave Dissipation; Na layer structure; Na 


flux; Na production rate 


